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We show that the topological Majorana quasiparticles hosted by semiconductor nanowires much longer than 
the superconducting coherence length are adiabatically connected with discrete zero-energy states generically 
occurring in short nanowires. We demonstrate that these zero-energy crossings can be tuned by an external 
magnetic field and are protected by the particle-hole symmetry. We study the evolution of the low-energy 
spectrum as a function of the magnetic field, wire length, and chemical potential, manifestly establishing that 
the low-energy physics of short wires is directly related to that occurring in long nanowires. In the presence of 
finite energy resolution, invariably operational under experimental conditions, the near-zero-energy states of a 
short wire give rise to a zero bias conductance peak over a finite magnetic field range. 



PACS numbers: 03.65.Yz 

The theoretical predictions ifTUTTl that proximity effect in- 
duced by ordinary s-wave superconductors, along with spin- 
orbit coupling and Zeeman spin splitting, could give rise to 
topological superconductivity have led to an intensive experi- 
mental search for Majorana fermions. Following specific the- 
oretical predictions [7, 8 ], a series of recent experimental pa- 
pers |[T2lfT5l have presented evidence for the existence of Ma- 
jorana modes in quasi- ID semiconductor nanowires. This ex- 
citement is further enhanced by the fact that these Majorana 
end modes can, in principle, be used to carry out fault-tolerant 
topological quantum computation (l6l[T7l, as envisioned orig- 
inally by Kitaev 1 18 ] more than 10 years ago. 

Any observation of the Majorana mode in solid state mate- 
rials is a rather important experimental discovery, therefore it 
is legitimate to ask critically whether the recent experimental 
findings are truly consistent with the theoretical predictions 
for the elusive Majorana particle. This is particularly impor- 
tant in view of the fact that the current experimental observa- 
tions (except Ref. fT3ll ) are based entirely on the existence 
of zero-bias conductance peaks (ZBCPs) in the differential 
tunneling measurements, which represents a necessary con- 
dition [ 9, 19 - 22 ] for the existence of the Majorana mode. The 
sufficient condition necessitates an interference experiment 
establishing the non-Abelian nature of these modes, which 
has not yet been performed. Since ZBCPs arise quite com- 
monly in both superconductors and semiconductors, it is of 
critical importance to carefully analyze the various experi- 
mental data to see whether the ZBCP is indeed consistent with 
the existence of the Majorana or is arising from other, pre- 
sumably more mundane, physical mechanisms. In addition, 
in short wires with lengths comparable to the superconduct- 
ing coherence length, it is commonly believed that the two 
end Majorana modes should hybridize and move away from 
zero bias |23] [24). Since disorder-induced chemical potential 
fluctuations could in principle break up a quantum wire into 
shorter segments (i.e., random charged puddles that would be- 
have as quantum dots), it is important to correctly identify the 



nature of the ZBCPs in small dots and to determine the fea- 
tures that distinguish them from Majorana ZBCPs. This is- 
sue has become urgent because the original observation of the 
ZBCP in long (> 2/im ) InSb nanowires by Mourik et al. lfT2l 
has recently been qualitatively reproduced in short (< 0.5/im) 
InAs nanowires by Das et al. fl5l , thus raising the important 
question of whether the ZBCPs in long and short wires are 
manifestations of the same qualitative physics or not. 

The goal of the current work is to critically investigate the 
wire length dependence of the ZBCP in superconducting (SC) 
nanowires and to clearly identify the nature of the ZBCP in 
short wires and its possible relationship to the Majorana zero 
modes in long wires. We establish that, for appropriate val- 
ues of the magnetic field, the lowest energy mode of the SC 
system is characterized by an adiabatic continuity as a func- 
tion of wire length and that ZBCPs generated by this near- 
zero-energy mode may exist even for wire lengths compara- 
ble to the SC coherence length. Therefore, although in short 
wires the whole notion of a topological phase with non-local 
zero-energy Majorana modes becomes meaningless, the zero- 
energy crossings associated with the ZBCPs can be viewed as 
signatures of remnant Majorana physics carried over from the 
long- wire topological phase. We thus believe that the ZBCPs 
observed in Refs. lfT2l and 031 for long and short wires, re- 
spectively, are adiabatically connected and, in some sense, are 
both manifestations of the predicted Majorana quasiparticles 
in topological quasi- ID superconductors 00. By solving 
numerically an effective tight-binding model for multiband 
semiconductor (SM) nanowires with realistic parameters, we 
calculate the energy spectrum and the density of states as func- 
tions of the wire length and externally tunable parameters - 
magnetic field and chemical potential. We show that robust 
zero-energy crossings generically occur in short nanowires 
at discrete values of the magnetic field B. When isolated 
in some parameter space, these accidental zero-energy cross- 
ings are protected by the particle-hole symmetry which is a 
hallmark of spinless p-wave superconductivity. With increas- 
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FIG. 1. (Color online) Low-energy semiconductor spectrum in the 
non- superconducting phase (orange/light gray) corresponding to the 
effective parameters of Eq. {2}. The black dots represent the dis- 
crete energy values for a short wire with L x & 400nm and periodic 
boundary conditions. 

ing wire length, the discrete values of B that characterize 
the short- wire zero-energy crossings are more dense and be- 
come continuously distributed over a certain range in the limit 
L — > oo, which corresponds to the topological SC phase that 
hosts the non-local Majorana zero-energy modes. We demon- 
strate that the zero-energy crossings in short wires generate 
ZBCPs that may look similar to those generated by Majorana 
zero-energy modes due to the limited experimental resolution. 
Our key finding is that particle-hole symmetry and the ex- 
ternal magnetic field together quite generically give rise to 
zero energy modes in ID spinless topological superconduc- 
tors, thus adiabatically connecting the ZBCP in a short wire 
and a long one. 

We consider a semiconductor (SM) nanowire with rectan- 
gular cross section L y x L z = 50nmx60nm and different 
length values along the wire L x . In the limit of an infinite 
wire, L x — » oo, the Hamiltonian describing the nanowire 
reads 



[/i - Y(j x - ak(jy}5 nrn - ia y q nrn cr x 



(1) 



where k = k x is the wave number, Oi are Pauli matrices as- 
sociated with the spin degree of freedom. In Eq. Q n = 
(n y ,n z ) and m = (m yj m z ) label different confinement- 
induced sub-bands described by the transverse wave func- 
tions 4> n (y) oc sm(n y 7ry/L y ) sm(n z 7Tz/L z ), e nm describes 
the SM spectrum without SO coupling, \i is the chemical 
potential, and T = g* (1bB/2 is the external Zeeman field 
along the x-direction. The effective parameters of the model, 
e.g., a, e nm , are calculated numerically following the proce- 
dure described in Ref. l25ll and using bare parameters for 
InAs (effective mass m e fj = 0.026mo and Rashba coeffi- 
cient aR = 0.2eVA). The low-energy spectrum of the wire 
is shown in Fig. [T] For a finite nanowire, the spectrum con- 
sists of discrete energy levels, see Fig.[T]where we have shown 
the spectrum for a short wire of length L x — 400nm. 

Next, we consider the SM nanowire proximity-coupled to 
an s-wave superconductor. The superconductor can be de- 
scribed by the BCS density of states v(E) = u F Q(\E\ - 



|A | 



where vp and A are the normal density of 



states at the Fermi level and the SC energy gap, respec- 
tively. By integrating out the SC degrees of freedom and 
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FIG. 2. (Color online) Low energy BdG spectrum as function of the 
chemical potential for finite wires of different lengths. The blue (dark 
gray) lines correspond to a vanishing Zeeman field, while the orange 
(light gray) line represents the lowest energy state for V = 0.18meV. 

linearizing the frequency dependence, one arrives at an ef- 
fective low-energy description of the system valid at ener- 
gies E <C A 1 25 ]. The corresponding BdG Hamiltonian for 
quasi- ID nanowire reads 



H nm (}v) — [fnm(^) fJ j 5 nrn \T z -\- YS nrn (J x T z 
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Here are Pauli matrices associated with the particle-hole 
degree of freedom. We have used the basis (u^ : u^v^^v^) 
for the Nambu spinors. Henceforth, we assume that the ef- 
fective SM-SC coupling is 7 = 75/ieV, which determines an 
induced SC gap Ai n d = 7Ao/(7 + Ao) = 50/ieV. For these 
parameters, the level spacing between different n x states be- 
comes larger than the SC gap A = 150/ieV in wires with 
L x < 0.5/im. In this regime, changing of the chemical poten- 
tial leads to significant variations of the energy corresponding 
to the lowest BdG state and of the number of quasiparticle 
states within the SC gap Ao. This behavior is illustrated in 
Fig. [2] In the absence of a Zeeman field (blue/dark gray lines), 
the minima of the BdG spectrum roughly correspond to the 
quantized energy levels E nxnyTlz = /1, with (n y , n z ) = (1, 1), 
i.e., the lowest energy band in Fig. [T] and different n x values. 
A similar behavior can be observed when the chemical poten- 
tial is in the vicinity of other band minima, e.g., \i — 18meV 
+A/i for the band with {n yi n z ) = (2,1), plus extra contri- 
butions from the lower energy bands. In the presence of a 
Zeeman field, the energy of the lowest-energy state decreases 
and eventually vanishes at a certain /i-dependent value of T. 
Note that, as a result of spin-orbit coupling, states with low n x 
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FIG. 3. (Color online) Low energy BdG spectrum as function of 
the Zeeman field for finite wires of different lengths. In the limit 
L x —> oo a Majorana zero mode appears above a critical field 
r c « O.lmeV. In finite wires, the mode acquires a finite energy due 
to the overlap of the states localized at ends of the wire. In very short 
wires (e.g., L x = 0.2/xm) the lowest energy state depends almost 
linearly on the Zeeman field. States characterized by different values 
of n x are coupled by the Rashba interaction and, consequently, the 
dependence of their energy on F is nonlinear. The chemical potential 
is ji = 18meV (bottom of the third band in Fig[TJ. 

depend strongly on the Zeeman field, while high n x states are 
weakly T-dependent. 

The dependence of the quasiparticle spectrum on the ap- 
plied magnetic field for several values of L x is shown in Fig [3] 
The lowest energy mode (red lines) is characterized by dis- 
crete zero-energy crossings that are robust against disorder, 
which we checked explicitly. In spinless superconductors, 
such isolated crossings are quite robust against perturbations 
due to the particle-hole symmetry, as we argue below. Con- 
sider k ■ p perturbation theory near a crossing point. The two 
zero-energy solutions and are related by particle-hole 
symmetry: = t x S&q. In order to open gap in the spectrum 
near the crossing point, the off-diagonal matrix element has to 
be non-zero: (^o|^|^i) ^ 0. Here V is a generic perturba- 
tion that satisfies particle-hole symmetry 

t x Vt x = —V T (3) 

However, using particle-hole symmetry we have Voi = 
(*o|^|*i> = /dztfgytfi = -fdx^r x V T ^ = 
— J dx^\V T ^Q = 0. Thus, particle-hole symmetry ensures 
the robustness of isolated zero-energy crossings. Another way 
of understanding the robustness of an isolated zero-energy 
crossing invokes fermion parity - one can show that the two 
zero-energy states and actually correspond to a dif- 
ferent fermion parity Q. However, the position of the zero- 
energy crossing point is non-universal and changes with the 
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FIG. 4. (Color online) Profiles of the lowest-energy states in 
nano wires of different lengths. The top panel shows a Majorana 
bound state localized near the end of a long wire (L x ^> £). In the 
other panels the red (dark gray) lines correspond to V = 0.05meV 
and the yellow (light gray) lines are for V = 0.18meV (see Fig. [3}. 
Increasing the Zeeman field mixes states with different values of n x 
and generates modes that tend to become localized near the end of 
the wire. This mechanism is absent in very short wires (bottom panel, 
(L x <C O) due to the wide energy separation between the quantized 
levels (see Fig' [2j. 

perturbation, since the diagonal matrix elements are non-zero 
(^ol^l^o) = -(^i|V|^i). In order to get rid of the zero- 
energy crossings one has to bring another pair of zero-energy 
states to the same point in the parameter space. Then, four 
states would hybridize with each other since two of them will 
now have the same fermion parity and eventually result in 
the avoided level crossings. This is illustrated in Fig. [3] for 
L x w 0.4/im. In this case, small variations of the chemi- 
cal potential will result in either two close zero-energy cross- 
ings (A/i < 0), or an avoided crossing (A/i > 0). These 
avoided level crossings would still, however, be near-zero- 
energy states, and may produce ZBCP in experimental sys- 
tems which invariably have finite energy resolutions. 

In wires with L x > £, the lowest-energy mode is sepa- 
rated from the excited states by a certain gap (see Fig. [3]). 
In the limit L x — oo, this gap becomes finite for T > T c , 
i.e., inside the topological SC phase, while the lowest energy 
mode becomes the Majorana zero mode. The corresponding 
states are localized near the ends of the wire and represent 
superpositions of different n x modes that are coupled by the 
spin-orbit coupling (see Fig. [4]). Note that, for T < T c the 
lowest energy state is delocalized. Remarkably, the tendency 
towards the formation of these low-energy localized states 
for appropriate magnetic fields is present even in relatively 
short wires with L x « £ (see Fig. [4]), which we interpret 



4 



> 

o 



> 



100 
50 


-50 
-100 
100 
50 


-50 
-100 
0.0 







Hi 










1 






r = 0.4meV| 


()OOoo°°°. 


i 








1 



0.5 1.0 



1.5 2.0 

L x (|um) 



2.5 3.0 3.5 



FIG. 5. (Color online) Evolution of the lowest-energy mode with the 
size of the wire for different values of the Zeeman splitting r = 0.2 
meV and r = 0.4 meV and chemical potential \i — 0. The shaded 
region corresponds to excited quasiparticle states. 

as remnant Majorana physics. However, when the spin-orbit 
coupling Esq = ^ e ff<^ 2 /2 is much smaller that the spac- 
ing between quantized n x levels, this mechanism becomes 
inaccessible (Fig. [4] lower panel). Thus, the adiabatic con- 
tinuity of the true Majorana-induced ZBCP in long wires to 
the near-zero-energy modes of the short wires requires strong 
spin-orbit coupling, finite magnetic field, and superconduct- 
ing pairing- exactly the same ingredients essential to produce 
the topological phase. 

The evolution of the low-energy spectrum with the wire 
length clearly illustrates the adiabatic continuity of the near- 
zero-energy mode, as shown in Fig. [5] Consider first the long- 
wire limit L x ^> £. Above the critical field T c > O.lmeV, 
the system is driven into a topological phase with Majo- 
rana zero-energy end states. In a finite system, the split- 
ting energy SE between Majorana modes has an oscillatory 
pre-factor, in addition to an exponentially-decaying envelope, 
SE oc smftpLx) exp(— L x /£) (23). The change of the sys- 
tem size or the magnetic field, which in turn changes kp and 
£, results in oscillations of the energy splitting, as shown in 
Fig. [5] and Fig. [3] With deceasing L x , the number of oscilla- 
tions within a given T interval decreases, while their ampli- 
tude increases, so that the gap separating the lowest-energy 
mode from the excited states collapses or SE exceeds Ao. At 
this point, all remnant features of localized Majorana modes 
completely disappear. 

Our results presented in Figs. [2|j5] clearly establish that 
robust near-zero-energy modes are generic in quasi- ID 
nano wires in the presence of spin-orbit coupling, magnetic 
field, and SC pairing, independently of the wire length. One 
of the possible experimental implications of these findings is 
illustrated in Fig. [6] In the presence of a finite energy reso- 
lution, the near-zero-energy mode is converted into a contin- 
uous ZBCP as a function of the magnetic field, as observed 
experimentally. This behavior is reflected by the dependence 
of the calculated density of states (DOS) on the applied mag- 
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FIG. 6. (Color online) Nearly zero-energy peak in the density of 
states (DOS) of a short wire with L x = 0.4/xm as function of the 
Zeeman field. The chemical potential is \i = 18meV. A decrease of 
the SC gap Ao with the increasing Zeeman field (by up to 60%) and 
a finite energy resolution of 10/xeV are included. 

netic field (see Fig. [6]). Note that the finite width of the zero 
energy peak in the DOS is determined by temperature and, 
more importantly, by the spin-orbit coupling of the near-zero 
mode to states outside the SC gap A . In short wires, where 
only a few states have energies inside the SC gap, this may be 
the dominant contribution. Coupling to metallic states from 
the leads could provide additional broadening. The apparent 
ZBCP will eventually split off and may come back again at 
still higher fields. Note that a Majorana zero mode could arise 
from coupling a finite number of quantum dots in a chain, as 
has recently been discussed in the literature l26l . 

We conclude by emphasizing our qualitative new finding 
with important implications for the current experiments prob- 
ing the existence of Majorana modes in hybrid semiconduc- 
tor structures. Except for very short wires characterized by 
quantized level spacings much larger than A and Esq, the 
system supports an adiabatically continuous low-energy mode 
that smoothly crosses over from a pseudo-Majorana regime 
with discrete zero-crossings to a zero-energy Majorana mode 
as the wire length increases. This mode is generically associ- 
ated with zero bias conductance peaks that, in a finite resolu- 
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tion measurement, extend over a finite magnetic field range. 
This indicates that the recent observations in Delft by Mourik 
et al. in long InSb nanowires and in Weizmann by Das et 
al. in short In As nanowires are adiabatically connected and 
are likely to be the expected signature for the predicted spin- 
less p-wave superconductivity characterized by the existence 
of the Majorana quasiparticles in long nanowires. 
This work is supported by Microsoft Q. 
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